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Abstract 

We calculate explicitly the space dependence of the radiative relaxation rates and associated level 
shifts for a dipole placed in the vicinity of the center of a spherical cavity with a large numerical 
aperture and a relatively low finesse. In particular, we give simple and useful analytic formulas for 
these quantities, that can be used with arbitrary mirrors transmissions. The vacuum field in the 
vicinity of the center of the cavity is actually equivalent to the one obtained in a microcavity, and 
this scheme allows one to predict significant cavity QED effects. 
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I. INTRODUCTION 



Many theoretical and experimental work has been devoted during recent years to the 
so-called "cavity QED" regime, where strong coupling is achieved between a few atoms 
and a field mode contained inside a microwave or optical cavity. In particular, it has been 
demonstrated that the spontaneous emission rate of an atom inside the cavity is different 
from its value in free space Jj, YA, , 0, [ll|, [3]. This effect can be discussed 

from several different approaches, and here it will be basically attributed to a change of 
the spectral density of the modes of the vacuum electromagnetic field, which is due to 
the cavity resonating structure. This approach is particularly convenient when the cavity 
does not have one single high-finesse mode, but rather many nearly degenerate modes, as 
it is the case in confocal or spherical cavities. More precisely, we will show that a "wide 
aperture" concentric resonator using spherical mirrors with a large numerical aperture, can 
in principle change significantly the spontaneous emission rate even with moderate finesse. 
Similar result were already demonstrated, using either a spherical cavity or "hour- 



glass" modes in a confocal cavity 



12 1 . In such experiments, the dipole has to sit within the 



active region volume, which is usually of very small size (of order (10A) 3 to (100A) 3 ). In 
refs [ill [l^ j a possible solution was implemented by reducing the cavity finesse in order to 
have an extended area in which a spherical wave is "self-imaged" on itself. However, getting 
large effects will put more severe constraints both on the quality of the cavity and on the 
localisation of the dipoles. 

A good understanding of these effects requires first to know the full space dependence 
of the cavity-induced damping and level shifts. In this paper, we will look at the situation 
where an atomic dipole lies close to the center of a spherical cavity with a large numerical 
aperture. We will show that large changes both in the atom damping rate and in its energy 
levels can be expected, even with a moderate cavity finesse, provided that the atom sits 
(relatively, but not extremely) close to the cavity center. In the following, we will assume 
that the cavity damping rate k is much larger than the free-space atom damping rate 7 mc - 
In that case, the cavity still acts as a continuum with respect to the atomic relaxation, and 
the damping rate and level shift of an atom at point r are given by [l^ |: 

W = S£(d.e fe (r)) 2 <5(u; fe -u, ) (1) 
n k 
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(2) 



where the summation are taken over a complete set of modes denoted by the index k. The 
resonance frequency and field at point r for mode k are respectively denoted Uk and efc(r), 
while the atom resonance frequency and dipole are respectively u and dD, where D is 
a dimentionless combination of raising and lowering atomic operators. We note that the 
free-space value of A(r) is a diverging quantity, which is usually assumed to be absorbed in 
the definition of the atomic levels; therefore, one considers here only the (finite) change of 
A(r) with respect to this free-space value, that will be denoted A'(r) : 



The purpose of this paper is to present an explicit calculation of A'(r) and T(r)/2, or 
equivalently of the modification of the (3-D) vacuum modes spectral density due to the 
presence of the cavity. For definitiveness, we will consider the case of an "open" spherical 
cavity, of radius R and of reflectivity and transmittivity coefficients r and t, with r 2 + t 2 = 1. 
The cavity can be "open" in the sense that it is made of two separate concentric mirrors 
which do not cover all 4-7T steradians. We will assume that kR ^> 1 (typically kR = 10 5 with 
k = u>/c), and a moderate cavity finesse (in the range 10-100). These parameters seems 
accessible from an experimental point of view, and we will show in the next sections that 
they allow one to get quite significant cavity-induced effects. 

II. MODES OF A LARGE APERTURE CONCENTRIC CAVITY 

We shall first consider the formal case of a scalar field, before turning to the real transverse 
electromagnetic field. Following the ideas of scattering theory, we propose a computation 
scheme where the propagation equations are cast in a form that is suitable for the determi- 
nation of the mode structure and that allows a convenient ray-optics formulation. 

A. Propagation of a scalar field 

We look for stationnary solutions (j)(r,t) = 4>(r)e~ lu)t of 



A'(r) = A C(ro (r) - A mc (r) 



(3) 



A0 - - 2 dU = 



(4) 
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that is, in spherical coordinates (r, 9, (j>) 



^9 r 2 r (r0) + A; 2 0+lA s = O (5) 



where we note k = u/c and introduce the spherical laplacian 



<9<m + + — (6) 



s sin 2 ^ ' ~ w ' sinO' 

with eigenvalues — 1(1 + 1), I > 0. 

In the far-field regime: r — > oo, we can write unambiguously 

0(r = rfl) = — /-'(A) + — /J^) (7) 
r r 

In the following we will often omit the argument Q of J out ' in . The quantity /° ut,m depend 
slowly on r: <9 r /° ut,m ~ /° ut ' in /r and tend to large r angular distributions: 



/r ,m — /r ,m (s) 

Such an 'out' field occurs for instance in the case of a radiating localized source: its squared 
amplitude then corresponds to the power radiated along Cl per unit solid angle. Here we 
also allow for incoming radiation, focused on a localized region — which, if not absorbed, 
turns after focusing into outgoing radiation. Note that this separation in 'in' and 'out' field 
is not possible too close to the origin, when the Poynting vector is no more almost radial. 

1. Far- field solution 

We have separate propagation equations for /° ut ' in : the 'in' field obeying 

drfr + ^f&fr + ^sfr = (9) 

Defining 5f — f r — we obtain 

d r Sf = ^sAsfoo ~ ^sSf - j- k 8l5f (10) 

with orders SJ - H ¥2 ¥2 

r kr A kr A kr z 

(we are only interested in the large r asymptotics, so we first ignore A s to get the orders of 
the different terms). 
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The first term in the r.h.s. gives the leading behaviour 

and a systematic expansion can be obtained by solving iteratively (fTUjl producing terms 

A? [ n > 1 

:/oo, " (12) 



(Ar)"+^^ p > 

In the sequel, we shall be interested in the field at ~ 100A off the origin: in geometric 
optics this involves light-rays with an impact parameter smaller than 100A, or photons with 
orbital angular momentum smaller than 100/L Therefore, in the multipole expansion of the 
field, we only keep harmonics with I < 100 : A s is now at most of the order 10 4 . Using our 
typical value kR ~ 10 5 , the terms f!12D with p > 1 are then negligible, and all p = terms 
are obtained by neglecting the last term in (jlOj) . which is equivalent to replacing @ with 

= ^sfr (13) 

Its solution 

fr=e iA e /2kr foo (M) 

then gives all the p = terms of the expansion ()12|). It is easy to show that with A s ~ 
10 4 , kr ~ 10 5 the magnitude of this first term is a few percent, while the next two terms 
range as 10~ 7 . Therefore, we shall use only the first term = in the following. 

To conclude this section, we extend these results to the case of outgoing waves: so far we 
only considered incoming radiation, but the analogue of (jHJ) is simply given by changing i to 
—i, and all results are easily transposed under complex conjugation, as an example of time 
reversal. In particular, we shall use 

/r = e-%r (is) 



2. Solution at any distance 

Up to now, we only studied asymptotic expansions of solutions of the wave equation in 
spherical coordinates, expressing f r (fl) for large r in terms of the values taken by foo in the 
neighbourhood of Cl. We will now derive an integral equation giving f r at any finite distance 
in terms of the function f^. 
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In the case of a stationary wave propagating in free space (with sources at infinity), we 
expect the knowledge of /™ to enable us to determine the solution <ft(r) everywhere, and in 
particular its asymptotic behaviour f™ tm . forming wave-packets, this amounts to constructing 
the solution at any time, given its value in the infinite past. We claim that the exact solution 
is obtained by a sum of plane waves 

0(r) = -lik J ^f«(£l)e^> (16) 

/™(fi) being the amplitude of the wave coming from direction Q with wave- vector — kCl. 

Obviously, the proposed solution does satisfy the wave equation, as a superposition of 
plane-waves; to prove our statement it thus suffices to verify that has the right asymp- 
totic behaviour /™. But for large r the integral is dominated by its points of stationary 
phase; suppose, for definiteness, that the axis 9 = is in the direction of r: then the phase 
in (|16|) — ikr cos 8 is stationary at 9 = and 9 = tt. Near each of these points, the lead- 
ing contribution to the integral will be of order 1/kr, with corrections corresponding to 
higher powers of 1/kr: so, the asymptotic part m ' out is entirely determined by the leading 
contribution to the integral. Using 



dQe- lAcosd ~ ^e" lA (17) 

near 9=0 A 



A^oo 



we obtain the contributions of neighbourhoods of 9 = and 9 = tt to </>(r) 



p —ikr p ikr 



f»(9 = 0) ; fZ{6 = tt) (18) 

r r 

corresponding respectively to the 'in' and 'out' fields, as could easily have been figured out. 

We recognize the right 'in'-field in this expansion, and have proved the validity of (fT7)|) . 
Moreover, we have obtained the following relation between 'in' and 'out' fields 

= -O-fi) (19) 

If a wave is focused, it emerges in the opposite direction, with the opposite phase. We also 
recognize in the i factor in the integral : 

0(r) = ^ f dhf»{h)e*- kCl > (20) 



the relative 7r/2 phase at the focus point. 

Corrections to this leading behaviour will produce the asymptotic expansion of /™' out 
in powers of 1/kr, the neighbourhood of 9 = contributing to /™ and that of 9 = tt to 
f° nt . In this way, one can recover, with longer calculations, the results of the preceding 
section. For instance, including the first correction to f* n amounts to replacing f™(9 = 0) 
withf£(0)+i/2kr (A s /£)(0). 



3. A complete set of explicit solutions 

We know that the angular distribution of the field has variations with r given by an opera- 
tor expressed with A s : thus, eigenf unctions of A s - spherical harmonics - give r-independent 
angular distributions (up to normalization and phase), that is, factorized solutions: 

Y lm {ti)j{r) (21) 

The wave equation and the smoothness of the solution at r = then determine j(r) up to 
a constant factor: 

3{r) = 7^— (22) 



Using the asymptotics of Bessel functions 

J l+1/2 {kr) [2 1 . , irl 1(1 + 1) 



kr y 7i kr \ 2 2kr 



we obtain the large r behaviour of this explicit solution 



^ - (* - i + ^) + (23) 



(f)(r = rn) = -2ik(-i)\l -Y lm (n) t+1/ ll ' (24) 

° ykr 



Y lm (Q) — -e-^ - (-l) l Y lm (Q)— -e^ + 0{l/r 2 

We can verify in this particular case the general relation f^iVL) = — Cl) and check the 
action of expiA s /2kr on /™ to the accuracy of (j24j) . 

We can also check the expression for the field at finite distance (J16|) : noting that /™ = 
and choosing normalized spherical harmonics 

dii , , „ 2 



(^|^m)= / 4^l y ^l =1 ( 25 ) 



we decompose 



and obtain 



fZ = E( Y lm I f£)Ylm (26) 
Im 

(i) = ew~ i a^f 2 yun) J -^l (27) 



The rotationally invariant quantity J2 m ^Zm(^)^m(^') is conveniently evaluated when the 
axis of reference is chosen along Cl and has value (21 + l)p(cosa), P being the /* h Legendre 
polynomial, and a the angle between fl and fV. Using the formula 

v / W/ ^ r) (2I + l)P(cosa) = ^e lkrcosa (28) 
i>o V^r V vr 

we finally obtain 

0(r) = ( W^' r | fZ(n'))n> (29) 

which reproduces (JTHJ). 



B. Modes of concentric cavities 

1. Perfect spherical resonator 

The explicit solutions given above allow us to determine the field modes inside a perfectly 
reflecting sphere, with radius E>i when we classify modes according to their spherical 
symmetry (quantum numbers l,m), the requirement that the field shall vanish on the inner 
face of the cavity 

f R = fRe- lkR + fTe lkR = (30) 
reads J l+ y (kR) = 0, or according to (J23|l : 

^-f' + ^T = ° t moc H ( 31 ) 

for the mode I, m, where we have written all significant terms for I ~ 100; kR ~ 10 5 , 
obtaining by the way the lowest order for which the /-degeneracy is disproved. The eigen- 
frequencies are then 

kc c ( I 1(1+ 1)\ , nn . 
n being the number of radial nodes. 
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With our numerical values the /-frequency shift has relative magnitude which is up to 
l/2nkR ~ 2.10" 4 , and is therefore very small compared to the cavity linewidth. We note 
that at a point close to the center, small / modes are more important, since photons travelling 
close to the origin carry a small orbital momentum. This point will be made quantitative 
later, when we will discuss the case of a spherical cavity with finite transmission. 

2. Modes in an open cavity 

a. Principle of the determination of the mode structure 
We recognize the vacuum fluctuations in a concentric cavity as induced by the vacuum 
fluctuations of the outer void space which enter into the cavity: we will thus start our 
analysis by studying how any incident radiation can enter in the open concentric resonator. 

We consider two spherical mirrors, facing each other in vacuum, with common center O 
- as if in the preceding example the tropical zone of the sphere were transparent, while the 
polar zones remained coated. For our computation, we shall replace the infinite vacuum 
with the inner volume of a very large sphere centered at O, thus replacing a true continuum 
with a discrete series of very closely spaced lines. We will use 1Z to denote the radius of the 
outer closed sphere, and R for the inner sphere, partially covered by mirrors; and use the 
following notation for the field of an eigenmode 

TZ>r> R : - e tkr f° ut + -e~ lkr f r n 
r r 

R>r : 0(r) (33) 

Note that the decomposition between 'in' and 'out' fields in the first equation is allowed by 
our choice to study only modes which contribute to <f> near O, that is, not too unfocused. 

We can formally extend for larger values of r as if there were no cavity at all, and write 
the far-field decomposition 

-^ kr gT + ^e~ ikr g™ (34) 

g™ is the incoming radiation which induces in void space the same field near O as /™ does in 
the presence of mirrors. We are ensured that g^ l (Q) = —g™(—£i) since that field propagates 
through the origin. However, the same equality does not hold for /: to understand the 
relation between f m and / out , we note that the incoming radiation f m induces a field in the 
open cavity, which in turn (perhaps after resonance) emits an outgoing field / out ; writing 
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the precise relation would require to solve the propagation equation for r > R, write the 
boundary conditions on the mirrors, and obtain the condition on f ul < ont for the existence of 
a solution between the mirrors satisfying g^fa) = — g™(— Cl). In some sense, g directly 
goes through the origin, while f m turns into / out after reflection on the mirrors, or multiple 
reflections in the cavity. We will not try to write any explicit formula for that, but use some 
general properties of the relation 

f» ^ fT = SfZ (35) 

in close analogy to scattering theory. 

We will assume that the losses on the mirrors are negligible when compared to the trans- 
mittivity: the balance between the incoming and outgoing energy fluxes 

\\fT\\ 2 = \\a\ 2 = J§\fz\ 2 06) 

requires S to be unitary. 

We now express the condition of perfect reflection on the inner face of the large closed 
cavity, i.e. that the field vanishes for r = 1Z : 

Once again, we are interested only in those modes that are focused enough so that they may 
contribute to the field near O: I being bounded, if 1Z is large enough we may use f n ~ 
and rewrite the preceding formula 

SJZ = ~e- 2lkn fZ (38) 

and formulate the problem of determining the modes of the open cavity (enclosed in a large 
one) as follows: find a wavenumber k such that — e ~ 2lkn is an eigenvalue of the unitary 
operator S k and identify the corresponding eigenvector. We precised our notation and used 
Sk to make clear the dependency of the scattering operator on the frequency at which the 
open cavity is excited. However, the problem is not so intricated since Sk depends slowly 
on k on the scale of the free spectral range of the 7£-cavity: 

5kRF<£l S k ~Sk + 6k (39) 

with T the finesse of the i?-cavity. So, in order to find the eigenmodes with wavenumbers 
k ~ k we have to find the eigenvalues and eigenvectors of S ko ; then adjust k so that —e~ 2%kn 
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coincides with any chosen eigenvalue of §k - Clearly, the mode structure so obtained will 
be (locally) periodic: the same eigenvector of Sk occuring as the far-field of a mode every 
5k = it/11 or Sv = c/2K. 

We would like to stress the close analogy between this general scheme and scattering 
theory. The operator S is to be thought of as an S-matrix in interaction representation: 
indeed, we can actually take the limit of large r for the connection between f™ and /° ut , but 
not for the relation between e~ tkr 'f" 1 and e lkr f° nt . Following our analogy, we may say that 
the radial propagation with phase factor e tkr corresponds to free evolution in perturbation 
theory, while the orthoradial propagation of light with changes in the angular distribution 
/ corresponds to the perturbation and asymptotically vanishes (at large r vs. at large times 
in scattering theory). 



b. Normalization of the modes 



We define the total energy of the field to be / d 3 r | 0(r) | 2 = 2 / d 3 r 2 eal im ° where rea i(r, t) 



Re 



. For large enough 1Z this energy integral is dominated by large r regions where 



we can approximate 

0(r) ~ ie^/r + \e- ikr fZ (40) 

and so obtain for the energy 

^ = 47r^(||^|| 2 +||/rH 2 ) = 87r^||^|| 2 (41) 

We decide to call vacuum the state in which every mode is excited with energy 1 (in real 
electrodynamics we should use Hcu/2) so that the normalization of any mode in vacuum is 

1 1 f oo 1 1 vac SllTZ- (^^) 

We can easily express the fluctuations of the ordinary vacuum field (here, ordinary means 
in infinite space) or rather their spectral density: in a range of frequencies Su we have 

V r 4ttV 



/ rf 3 k= ^L„2 5u (43) 



(2tt) 3 hv 

modes for the case of a large volume V with periodic boundary conditions, while each mode 
contributes 1/V to | 0(r) | 2 since the energy of any mode is uniformly distributed in the 
volume. Finallyone has the expected result in ordinary vacuum: 

An 



\<K*)\% = ^6v (44) 
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In the next part, we will derive an analogous formula for the case of an open concentric 
cavity and compute by how much the vacuum fluctuations (near the center) are amplified 
or reduced by the presence of mirrors at about 1cm. Before that, a last comment is in 
order: to obtain the above expression for the ordinary vacuum field we used the fact that 
the contributions of different modes add up incoherently. This is always true when we use a 
basis of modes in which the energy operator (hamiltonian of the field) is diagonal: in general, 
eigenmodes are non-degenerate in frequency and this condition is automatically satisfied. 
This vanishing of the off-diagonal matrix elements of the energy reads 

|rf 3 r^(r)0 2 (r)=O (45) 

and expresses the orthogonality of different modes with respect to volume integral. 

The discussion of the preceding part enables us to assert a more precise property, which 
will be essential in the following: even the angular overlap of different modes is 0. Indeed, 
we argued that eigenmodes were obtained by diagonalizing the unitary operator Sk - but 
evidently, two different eigenfunctions of the same unitary operator are orthogonal. Thus, 
the far-fields of different modes have a vanishing angular overlap, unless the two particular 
modes do have the same far-field asymptotics : in the latter case, their number of radial 
nodes being different ensures the vanishing of their volume overlap. 

That property can be stated differently: if we consider all eigenmodes in a range of 
frequencies 5u = c/2H and associate to every such mode its /™ we obtain a complete 
orthogonal set of normalized angular functions. Should we consider a larger range 5u every 
member of this orthogonal family would then be counted ITLbv jc times. 

c. Modification of the vacuum field in an open cavity 
We will find an expression for the field at points close to the origin (r < a few hundreds A that 
is a few dozens pm). All light rays that pass so close to the center will then reflect almost 
normally on the mirrors: we will assume that the wave-fronts of all modes are sufficiently 
tangent to the mirrors for us to use the % = reflectivity and transmittivity coefficients p, r 
(p = 1 for a perfect mirror). 

Let us first compute the field induced in the open cavity by incident radiation /™: we 
will note f™ = U r f™ with U r ~ e * A s/ 2fer according to our previous results. A similar relation 
holds for 'out' fields with U r = or U~ x since U r is (almost) unitary. On the outer face 
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of the mirror, the field has incoming and outgoing amplitudes 

e ~zkR „j n 
— R— U RJoc 



(46) 



while on the inner face we find 



U R g 



m 



with 



R ^-Ki/oo 

(47) 

ikR 



stfW = -s£(-n) : g/Z? = -Pg& (48) 

P being the parity operator, that commutes with A s and U R . The g m wave has two contri- 
butions: partial transmission of f m and partial reflection of g ont : 

e~ tkR U R <£(fi) = -p(Cl)e ikR U R g£(Cl) + r(Cl) e - ikR U R f£(£l) (49) 

where we explicitly write the angular dependency of p, t: in particular, out of the mirrors 
p = and r = 1. We can easily solve ()48I49|) to obtain 

which yields a formula for the field at any point through 

0(r) = (2ike ikClr | «£(0)) fl (51) 

Introducing the shorter notation 

T = r{(l) ^U R (52) 

U 2 R -e^ kR Pp{Q) V ' 

we obtain for the vacuum field 

(i 0(r) | 2 ) vac = £ ( 2 ^ r I ^ + ^ I /£> </£ I ^ ^ I 2 ^ r > ( 53 ) 

modes: /jg 

The above-mentioned property of orthogonality of the modes now gives a considerable sim- 
plification since we do not need to know the precise expression of all the modes, but only 
their total contribution to the physically meaningful quantity (| 0(r) | 2 ) vac in a given fre- 
quency range: so, no matter what the true f™ may look like, they surely give a closure 
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relation. Recalling that the /™ associated to the modes in a frequency range of c/2TZ form 
a complete orthogonal set normalized according to (}4*2j) we see that 



modes in 5v. /™ 

and is actually independent of the large cavity we used to mimic the infinite vacuum. We 
then have 

(l0(r)| 2 ) vac = ^||^T.e^ r || 2 (55) 

7TC 

Using the unitarity of Ur and the expression of ordinary vacuum fluctuations we find 

(| 0(r) | 2 ) cav _ ,, -ifcQri|2 _ II \ iA a /2kR -ifcfir||2 (r^\ 

(I0(r)l 2 )vac ~ ~ l] e iA °/ kR -e 2ikR Pp " 11 { ' 

3. Numerical study and ray optics interpretation 
a. Numerical results 

As expressed in (J56|) . we shall evaluate the action of an operator on the function Q i— > e~ tknr 
and then compute the squared norm of the resulting function: this can be done numerically, 
decomposing functions on spherical harmonics 



e lkQr = E H)Vi 'T - nm(r/r)y,, m (fl) (57) 

Z,m;0<|m|<' ^ 

and computing matrix elements of the fc-dependent operator in the basis of spherical har- 
monics. 

We only considered axially-symmetric cavities: in that case, the operators involved con- 
serve m and the result is expressed as a sum of contributions from the different m sectors. 
Moreover, m = gives the single contribution to the field on the axis. In the latter case, we 
could use truncated systems of up to Z max = 300 spherical harmonics and study the effect 
of truncation: the result was constant within 1% for Z max > 100. In the general case (field 
fluctuations at points away from the symmetry axis of the cavity), we used Z max = 100 with 
any m to study the spatial and spectral dependency of vacuum fluctuations. 

The graphs shown on fig. Qwere obtained with concentric mirrors of uniform reflectivity 
p = 0.98 (giving an intensity transmittivity T ~ 4%) covering 30% of the Att steradian; this 
is obtained for an half-aperture angle 9 m = 45° for each mirror. The frequency is set at 
the resonance value at the cavity center. We note the rapid oscillations of the amplitude of 
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FIG. 1: Spatial variation of vacuum fluctuations on the cavity axis, obtained from a full numerical 
calculation. The horizontal axis unit is 1/k, and the vertical axis unit is the standard vacuum 
level. The amplitude reflection coefficient of the mirrors is taken to be p = 0.98, and the numerical 
aperture of the cavity is 0.7. The right hand side of the figure is a zoom close to the cavity center. 

vacuum fluctuations near the center of the cavity; right at the center, the obtained value at 
resonance agrees approximately with the usual rough estimate 4/Tx AQ miTT / 4n . However, at 
a few /im away from the center the enhancement effect is halved, and then decreases further 
on a scale of ~ 15A. The next sections will give support to a qualitative and quantitative 
formulation of these facts in terms of ray-optics. 



b. Case of a closed cavity: ray optics interpretation 
We come back to the case of a closed spherical resonator, now allowing a non-zero trans- 
mittivity of the mirrors (but still negligible losses, as we always suppose in this article). We 
can apply to that particular case the formalism we developed for open cavity, and recover 
the known modes involving spherical harmonics: 

1 JA s /2kRy _ * ~U(l+l)/2kRy , R g\ 



e iA s /kR _ e 2ikRpp l ' m e -U{l+l) / kR _ (— ]V p e 2ikR 

so that the vacuum fluctuations (normalized to 1 for usual vacuum) read 

J l+ i /2 {krf 

E | e -u(i+i)/kR _ (_i)Jpe2ifcfl |2 2 ( 2/ + V kr ^ ^ 

where we recognize the usual resonance factor, including the slight non-degeneracy of modes 
(increasing I: lower resonant frequency) and see the explicit spatial-dependency of /-modes. 
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Averaging vacuum fluctuations in a frequency range much larger than the 'free' spectral 
interval of the cavity we can replace the frequency- dependent factor by its average value: 1, 
and use the following sum rule 

ljM)^ = ' (60) 

We actually recover usual vacuum if the cavity is large enough (with free spectral range 
smaller than the experimentally used frequency bandwidth). 

If we are interested in the behaviour of vacuum fluctuations in the vicinity of the center 
(a few microns) where only small I modes contribute, we may assume that all modes are 
degenerate and use a common resonant factor: 

T 



-l) l pe 



2ikR 12 



(61) 



The sum rules for Bessel functions 



t:odd/cvcn 



then allow us to write the vacuum field as 



T (\ sin 2kr\ T (\ sin 2 Ax \ 

+ -T1Z- + m ■ U--T1Z- (63) 



|1 - pe 2ikR \ 2 \2 4Ax J 1 1 + pe 2ikR \ 2 \2 Akr 
At any point, we have two series of resonant lines, in which the vacuum noise is distributed 
with weights y 2 ± sin 2Ax/4Ax. At a frequency which is resonant for the center of the cavity, 
the spatial dependency of the vacuum field shows a reduction by a factor 2 when one moves 
away from the center, as was noted above in the case of an open cavity; we interpret this as 
the distribution of vacuum fluctuations on the two series of lines: right at the center only 
the / = mode appear, but the odd I modes have no other common nodes and share ~50% 
of the vacuum noise away from the center. The same conclusions can be formulated in terms 
of light-rays: noting that 

1 sin 2kr r d(l 9 , A 1 sin 2kr r dQ . 9 , A ,„ JN 
- + — = / — cos 2 kttr, = / — sin 2 kttr 64 

2 Akr J 4tt 2 Akr J Air v 1 

we may express the vacuum fluctuations (still neglecting the /-dependency of resonant fre- 
quencies) as 

||yj2 CaV = / ? (h ^^cos 2 (Mir) + Z_^ s in 2 (^.r)'l (65) 

(I0( r )l )vac ' An \\1 - pe 2tkR \ 2 v 1 \l+pe 2lkR \ 2 J 
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and propose the following interpretation: through the center of the cavity we may draw a 
ray in any direction; that ray reflects on the inner face of the cavity back onto itself and gives 
rise to a system of stationnary waves. For any such ray, the field oscillations (forced by the 
outside vacuum) may have maximum amplitude or a node at the center: correspondingly, the 
stationnary wave will have squared amplitude cos 2 kClr or sin 2 k&r at the point of interest. 

As for the contribution of rays that support a mode having an antinode at the cavity 
center, they have different phases away from the origin according to their direction and 
thus contribute in the average with weight l / 2 '- the former considerations on the positions 
of nodes of Bessel functions are now reformulated as positions of nodes of stationary waves 
along rays with different directions. 
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FIG. 2: Spatial variation of vacuum fluctuations on the cavity axis, obtained from a simple ray- 
optics calculation. The horizontal axis unit is 1/fc, and the vertical axis unit is the standard vacuum 
level. The amplitude reflection coefficient of the mirrors is taken to be p = 0.98, and the numerical 
aperture of the cavity is 0.7. The right hand side of the figure is a zoom close to the cavity center, 
which shows also a comparison with the full calculation (fig. ^). 



c. Spherical aberrations 
As it can be seen on fig. |21 the simple ray-optics analysis described above agrees approxi- 
mately with the results of the operator-based numerical computation in the vicinity of the 
center, but fails to describe the vacuum fluctuations away from that area. 

Indeed, in order to handle field properties at a point located at a distance r form the 
origin, we shall not consider anymore the effect of rays going through the origin, but rather 
that of rays going through this point : such rays miss the origin by a distance d < r, and 
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thus carry t < kr orbital momentum in % units. These rays do not close after one round 
trip : if we do not move too far from the origin, we may still assume that any point is 
imaged onto its symmetric after one reflection, and so onto itself after two reflections, but 
the (twice) reflected way is tilted by an angle 89 = ^ sin(0). However, as we are considering 
finesse values in the range 10-100, and r, R values respectively smaller than 10 2 and 10 5 A, 
we may reasonably neglect this tilt and and associate to any ray an average value 9. Note 
that 9 corresponds to i via i = krsin(9). since i is conserved at reflection on the mirror, 
due to the symmetry with respect to a radius, the change in the direction of the light ray is 
actually accompanied by a lack of re-imaging of the point back to itself. So, our assumption 
really consists in neglecting both the tilt and the lack of re-imaging, and correspondingly in 
keeping the spatial modulation cos 2 kr and sin 2 kr in the expressions given above. 

What cannot be neglected, however, is the relative phase with which the reflected light 
comes back to the initial point : for off-center rays, a round trip involves propagation on 
a distance AR + sin 2 (9). Consequently, the frequency-dependant factor acquire an extra 
phase term and becomes 2ikR + sin 2 (6 ) ). As it can be seen on fig. El this phase term is 
basically responsible for the decrease of the spectral density as a function of the distance. 
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FIG. 3: Same as fig. [2 now including spherical aberrations, so that the spectral density is now 
decreasing as a function of the distance. The result of the full calculation (fig. is shown for 
comparison. 



d. Boundary effects and diffraction losses 
The formulation in terms of light rays described above can be straightforwardly extended 
from a closed to an open cavity : to any ray we associate an enhancement factor, which is 
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frequency-dependant if the ray meets the mirror, and is unity if the ray misses the mirror, 
as well as a spatially-dependant term describing the intensity modulation of the stationnary 
waves. 

When the two mirrors are identical, we use the results known for symmetrical Fabry- Perot 
resonators and obtain : 

7T^H = /?( T — ft cos 2 (m.r) + T , kr2 — sin 2 (MIr)) 

(66) 

Should the two mirrors be different, the formula would be easily modified, just as in the case 
of a usual Fabry-Perot resonator (see Appendix A). In particular, the maximal enhancement 
at the center of a symmetrical cavity subtending a total solid angle f2 m is : 

(l0(O)| 2 )cav = ^oc T fl m 

(l0(O)| 2 ) vac " 4tt + |l-p|M7r 1 ) 

It is worth noting that the ray-formula without spherical aberrations given above corre- 
sponds exactly to the result of the more rigourous analysis, when one neglects As, i.e., if one 
takes e % <=« equal to unity. Though the main effect of As is accounted for by spherical aberra- 
tions, a small discrepancy remains : for p m = 0.98 and = 0.3, the numerical computation 
in the basis of spherical harmonics yields an enhancement factor of 29.2 at the cavity center 
and at resonance, while the ray computation gives 30.4 in the same conditions. We shall 
explain this small difference by diffraction losses : those rays that would be reflected near 
the edge of the mirror are actually lost due to diffraction and fail to do as many round-trips 
as the other rays. This second effect of As can be estimated by looking for an approximate 

.Ac 

inverse of the operator (e l w —p{9)), valid near the mirror edge. The result of this procedure 
is that one can still use the previous formula for any detuning and at any point, provided 
that the boundary value 9 m is decreased to 9 e ff — 9 m — 59, with 69 = 1/y kRT for symmet- 



rical mirrors, and 59 = l/JkR(l — p? v ) for non-symmetrical mirrors, p av being the average 
reflectivity (pi + p 2 )/2 of the two mirrors. Applying this procedure to the above example, 
we shall substract 

1 sin(fl m )( M T |2 - 1) = 1.2, (68) 



VkRT v " vv |l-p 
which is quite satisfactory since 30.4 — 1.2 = 29.2. 

The comparison of the ray calculation and of the complete one for an open cavity is shown 

on fig. HJ with the same parameters as for fig. As it can be seen, the agreement is very 
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good, and justifies a posteriori the assumptions which have been made. It can therefore 
be concluded that the main correction to the naive calculation is the phase error due to 
spherical aberrations, with some small correction from the edge diffraction losses. These 
corrections are enough to get the right answer in the conditions that we are considering 
(R ~ 10 5 A, r smaller than 100A). 
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FIG. 4: Same as fig. [21 now including both spherical aberrations and diffraction losses. The result 
is now in very good agreement with the full calculation (fig. which is also shown for comparison. 



e. Focusing defects 

The numerical scheme of the previous sections allows us to systematically study the effects 
of mechanical defects in the cavity, e.g. defocusing, or non perfectly spherical mirrors. For 
instance, one can reproduce an axial mispositioning of the mirrors by a length S, by adding 
an imaginary part to the mirror reflectivity : p(9) = p Qe 2lk ^/ R ' 2 + s ' 2 +' 2RS cos e-R _ p Qe 2ifc<5 cos e ^ 
The first effect of such a mispositioning is to shift the resonance frequency. Correcting for 
this shift, the second effect is to decrease the enhancement effect, which is typically halved 
for 5 = 40nm(k5 = 0.3) with the previous parameters. Here again, it can be seen that 
the ray formula gives the right answer. As it was discussed above for positions outside 
the cavity center, the main feature is indeed the phase shift after one reflection, which is 
correctly described by the modified value of p(9), while the tilting and non- imaging effect 
can be neglected. 
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C. Polarization effects 



The above considerations, which were done for a scalar field, can be straightforwardly 
extended to the case of a vector field, which is needed to describe polarization effects. Here 
we will skip the explicit operatorial formulas, and give only the results obtained in the ray 
optics approximation. As previously, this approximate solution was checked by comparison 
with the complete numerical calculation, and found to be in complete agreement with it. 
For a transverse vector field 0(r) and for two polarization directions ei and e 2 , the results 
obtained in the scalar case (eq. lo^j) are then changed into : 

((ei.0(r))* x (e 2 .0(r))) cav f dCl 3 



<( ei .0(r))* x ( ei .0(r))> 



vac 



/^|( ei . e2 -(a ei )(a6 2 ))x 



(^b^ cos2(Hlr) + ll + ^p sin2( ^- r) ) (69) 

From this equation, the cavity induced damping and level shifts can be obtained using 
eqs. (1) and (2) by integration over the frequency, which is straightforward for the damping, 
and requires contour integration for the level shifts (see Appendix B). Finally, the effect of 
the cavity can be described to a very good approximation by the following formulas : 



T mc i^47r 2 V v d ' ' \\1 - pe 2l( t>\ 2 V ; 1 1 + pe 2l( t> \ 



where the notation Cl describes a direction in space, while is a cavity detuning parameter 
that will be detailed below. As previously, these expressions have a straighforward inter- 
pretation, because they appear basically as integrals over the direction of light rays : in the 
integral over the directions, p is the mirror reflectivity for rays subtended by the cavity, and 
is zero for rays outside the cavity solid angle. The different factors appearing in the integrals 
are detailed below. 

The first factor under the integral corresponds to polarisation effects, taking into account 
the transverse character of the field. 

The second (resonance) factor is of the usual Fabry-Perot form, where <fi is the cavity 
phase shift which includes first a term O — u$R/c. As it was shown before, in order to 
obtain a correct result outside the cavity center, must include also a contribution from 
spherical aberrations, that is : (p = (fio + fc(r ~ 2 ^' r ' > ^ This second term corresponds to the 
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extra phase shift experienced by rays going through point r while propagating along the 
Cl direction. The resonance factor has obviously different expressions for the damping and 
the lamb shift terms, which correspond respectively to the active and reactive parts of the 
coupling. This is clearly apparent from the integrals of eqs.l and 2, which involve either 
a delta function or a principal part. In the first case, the integration is trivial, and yields 
the resonance term of eq. [7U1 while in the second case the result is obtained by contour 
integration, and gives the "dispersive" second term of eq. O 

The third term under the integrals is the stationnary wave pattern corresponding either 
to odd modes (which have an anti-node in the center and a cos 2 (/cf2.r) space dependence) 
or to even modes (which have a node in the center and a sin 2 (/cf2.r) space dependence). 

Finally, the integration over the mirrors is conveniently performed in spherical coordi- 
nates, by taking the z axis along the cavity axis, and varying the azimuthal angle 6 from 
to Q m irrar — @m- Improved accuracy (better than 1%) is obtained if one takes into account 
the fact that the rays which would be reflected near the edge of the mirror are actually lost 
due to diffraction and fail to do as many round-trips as the other ones. As before, this effect 
can be taken into account very simply by decreasing 9 m to 9 e ff = 9 m — 56, with 56 = ^j^ RT 
for symmetrical mirrors. 

The first results which can be obtained from the previous formulas are obviously the shift 
and damping at the cavity center, as a function of the atom-cavity detuning. For a dipole 
orientation parallel to the cavity axis, we obtain straightforwardly : 

T par (0) _ AQ vac n | sin 2 6> m A | AQ cav ^ cos6> m (l + cosfl^ T 
T mc " 4vr U+ 2 )+ 4tt 11 2 } \l-pe»*°\* 

A par(°) _ AQ cav _ cos6> m (l + cos6> m ) p sin(20 o ) 
T vac 4tt 1 2 )\i- pe 2i4>o\2 W 

while for a dipole orientation perpendicular to the cavity axis, we have : 

T perp (0) _ AQ vac sin 2 6> m AQ cav cos6> m (l + cosfl m ) T 
T vac 4tt 1 4 j 4tt 1 4 j 1 1 - pe 2i ^ p 

A per P (Q) _ AQ caV n | cos6> m (l + cos9 m ) ^ p sin(20 o ) 
We note that these expressions yield for a randomly oriented dipole : 

r a „(o) = Aa uac An cav t a'(o) = An cav P sin(20 o ) 

T mc 4vr 4tt |l-pe 2i ^|2' r mc 4vr \l-pe 2i ^\2 ^ ' 
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corresponding to the scalar case already given above. We note that these results are the 
same as those given in ref. [jj, up to factor two resulting from the fact that this reference 
was considering spatially averaged values rather than the peak value at the cavity center 
(see below for the space dependence). These functions are plotted on fig. EJfor = 0.3 
and p = 0.98. It can be seen that very significant effects occur for these quite reasonable 
parameters, yielding more than 30-fold increase in the damping rate at the cavity center. 
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FIG. 5: Normalized damping T(0) /T vac (left) and level shift A'(0) /T vac (right) at the cavity center, 
as a function of the atom-cavity detuning normalized to the cavity linewidth. The amplitude 
reflection coefficient of the mirrors is taken to be p = 0.98, and the numerical aperture of the 
cavity is 0.7. The upper curves correspond to a dipole oriented perpendicular to the cavity axis, 
and the lower curves to a dipole oriented along the cavity axis. 

The above formulas also give the damping and level shifts as a function of space for a 
given frequency, which are an important result of the present paper. The results in the 
most general case where the two mirrors have different reflectivities are given in Appendix 
C. Two atom-cavity detunings are specially worth looking at : the resonant frequency at 
the cavity center, which yields maximum change in the damping rate but no cavity shift, 
and frequencies detuned by plus or minus half a cavity linewidth, which yield maximum 

n 

cavity shifts. These results will be exploited in the following paper |l3j, which deals with 
vacuum-induced light forces acting on an atom close to the cavity center. 
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III. CONCLUSION 



As a conclusion, we have calculated explicitly the cavity induced damping and level 
shifts for an atomic dipole close to the center of a spherical cavity. Our results are valid 
for arbitrary (large) aperture and (not too large) mirrors reflectivities (for the most general 
case see Appendix C). These results show that macroscopic cavities with large numerical 
apertures are interesting candidates for cavity QED experiments in the optical domain. In 
particular, we show in a joint paper k| that the cavity-induced level shifts are responsible 
for a "vacuum-field" force on at atom moving close to the cavity center [14, 15]. An explicit 
expression of the trapping potential can be obtained from the results given above. 



Appendix A 

When the mirrors'transmission are different, eq. IHElcan be generalized to : 



(|0(r)P> w J 4?r > 2R 
where r = I r I, and : 



n 2 l + p 2 2 + 2p 2 cos(2(0- Mir)) 

M = } — — 

2 (i_ e -4 4 pip2 ) (i-<Ai* pxP2 ) 

t 2 2 (l + pi 2 + 2 px cos(2 (0 + Mir)) 



(78) 



2 (l-e- 4 ^p lP2 ) (l-e 4 ^ Pl p 2 ) 

where for each mirror r 2 + p 2 = 1. This equation can also be written in the less compact 
but more transparent form : 

M = (1 ^ lP2) ( l.j PlP2 + (Pl + / : 2) COS | 20)) cos 2 (AAr) + 



4 ^PiP 2 ) (l-e 4 ^p lP2 ) 

sin 2 (Mir) + 

sin(2Mlr) (79) 



'1 ~ Pi P2) (1 + Pi P2 - (pi + Pg) cos(2 )) 
(l-e- 4 ^p!p 2 ) (l-e 4 ^piP 2 ) 
(1 + Pi p-i) (P2 - Pi) sin(2 4>) 



(1 -e- 4 ^p!p 2 ) (l-e 4i *pip 2 ) 

which has the same interpretation as eq. EH1 : the sin 2 (A;f2.r) and cos 2 (fcflr) correspond 
to the contributions of the in-phase and out-of-phase standing waves along the direction 
Cl, while sin(2A;f2.r) = 2 sin(fcfi.r) cos(kQ.r) is an interference term due to the intensity 
inbalance between the forward and backward contributions. 
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In the case of a high finesse cavity (ti, t% « 1), these equations can be rewritten : 

2 (if cos 2 (0 - fcfi.r) + r| cos 2 (0 + Mir)) 
M = (e 4 ^ - 1 + (r 2 + r|)/2) (e" 4 ^ - 1 + (r 2 + r|)/2) (80) 

or alternatively : 

2 ( r i 2 + t%) (cos 2 (0) cos 2 (A;fi.r) + sin 2 (0) sin 2 (A;fi.r)) + (r 2 - rf) sin(2 0) sin(2 jfefij) 



M 



( e «* _ x + ( r 2 + r 2 )/2 ) ( e -4^ _ l + ( r 2 + r |)/ 2 ) 



(81) 



For a symetrical high-finesse cavity with T\ = T2 = r, one obtains finally 

4r 2 (cos 2 (</>) cos 2 (Mir) + sin 2 (0) s\n 2 (k(l.r 

M = i — — -L (82) 

(e 4 ^ - 1 + r 2 ) (e~ 4 ^ - 1 + r 2N 



which can also be obtained directly from eq. 



Appendix B 



Let us consider the normalized Airy function : 



^ ~ 1 ■ Fsin 2 " |l-pe 2 ^' 2 (83) 



where F is related to the mirrors amplitude reflectivity by F = 4p/(l — p) 2 . In order to 
calculate the level shift, we need to evaluate the principal part integral : 

A(0)= (84) 

Replacing C(4> — 5) by its uneven part \{C{4> — S) — £(<f) + 5)), A(0) can be expressed as a 
standard integral : 



AM _1 f dS VTTf VTTf 

A(0) " 2 J T ( l + Fsin 2 (5-0) " l + Fsin 2 (5 + 0) j (85) 

This quantity can be evaluated by contour integration, using the zeros of the denominators 
1 + F sin 2 (5 ± (f>), which are respectively 8~ = —(f> ± i(3 + nir, and 6^ = (f> ± i/3 + nir, where 
(3 > and sinh 2 (3 = 1/F. Using a contour in the lower part of the complex plane, which 
includes the poles 5„ = ±(f> — i/3 + nrr with n — ... — 1, 0, 1, we obtain for instance : 

rdS 1 _ -2m 1 _ 2tt 1 

J Tl + Fsin 2 (5-0) ~ Fsin(-2i/3) n <f> - i(3 + nir ~ Fsinh(2/3) tan(0 - i/3) ^ 
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Using sinh(2/?) = 2VTTF and tan(0 - iff) = agg^g^ , we obtain thus : 



AU] = 7T / 1 1 \ = vrF / sin(20) \ 

W 2Vtan(0-i / S) tan(-0 - iff) ) 2 1^1 + F sin 2 J 1 J 

Coming back to mirrors reflectivities, we obtain finally the expression used in eq. O : 



Appendix C 



In the case of asymetrical mirrors with amplitude transmitivities p\ and p 2 , one can use 
the formulas of Appendix B, changing F into F' = 4pip 2 /(l — P1P2) 2 , and into 0' = 20, 
so that : 

yv ! 2 \l + F' sin 2 07 v ; 

In addition, we need to evaluate the principal parts for £(0—5) multiplied either by cos(0— 5) 
or sin(0 — 5). Taking for instance the cosine part, we obtain : 

A c (0) = J £(0 - 5) cos(0 - 5) (90) 

This can be done as before, and we have : 

dS cos(<5 — 0) — 2in cos(— if3 + wr) 27r cosh(/9) 1 

T 1 + F sin 2 (5 - 0) = F sm(-2iff) n <p - i(3 + mr = Fsinh(2/5) sin(0 - iff) ^ ' 



Using sin(0 — iff) = (a/1 + F sin — i cos 0) / \/F we get : 



7T + F) 1 1 sin0 

Applying the same method for the sine part, we obtain finally : 



» / ,t\ sm0 . . ,,. cos 14) 

A c 0' =tt 1 + F y 2 A s 0' = -7T — -— 93 

1 + b ' sin 1 + b ' sin ^ 



From these formulas, we obtain the damping and level shift in the asymetrical case : 

T mc J 4* An 2 l l rf ;A (l-e- 4 ^p lP2 ) (l-e 4 ^p lP2 ) 1 ; 

(l-pip 2 ) (1 + P1P2 (pi + P2) cos(2 0)) gin 2^Q ^ _|_ 



(l_ e -4^ pipa ) ( 1 _ e 4 l pip2 
[1 + P1P2) (P2 - Pi) sin(2 0) 



(1 _ e -4i0 pip2 ) (l-e 4 ^pip 2 ) 



sin(2Mlr)} (94) 
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A'(r) 

r 

1 vac 



f dtt 3 ,d.£l^ ( Pl p 2 sin(40) + (pi +p 2 )(l + Pip 2 ) sin(20)/2 2 ft ^ 
y-br Itt 2 (1 ~ ( ~ )M (i-e-4^p lP2 ) (i_ e 4^ pip2 ) cos(/d2.r) + 

sin 2 (/cfir) + 



Pip 2 sin(40) - pi +p 2 )(l + PiP2) sin (20)/2 _ 2f 

(1 -e- 4 ^p!p 2 ) (l-e 4 ^p!p 2 ) 
(1 - pi p 2 ) (Pi - P2) cos(2 0)/2 



sin(2Ml..r)}. (95) 



(1 -e- 4 ^p!p 2 ) (1 -e 4 ^pxp 2 ) 
An particularly interesting case is a one-mirror cavity (pi = p, p 2 = 0), for which 

T(r) /■ dtt 3 dA 



r 

- 1 - vac 



A'(r) f dtl 3,„ ,d.Q.~. r p . > ,„„s 



r 



vac 



Neglecting spherical aberrations, one has as before <fi = ujqR/c, while r corresponds to 
the atom's position with respect to the mirror's center of curvature. We note that these 
equations have the correct behaviour T(r) = T vac and A'(r) = if p = (no cavity). In 
the case of a small solid angle e = Q/(4tt) subtended by the spherical mirror and a dipole 
orthogonal to the "cavity" axis Oz, one gets : 

« 1 + ^ cos(2(A;z + 0)), ^ « ^ sin(2(^ + 0)). (98) 

J- tiac ^ J- vac ^ 

From eq. UHl the total phase (kz + <p) corresponds to the distance I between the atom and 
mirror 1. Taking into account that kz + <p = 0, or equivalently kl = tt/2 mod 7T, are antinodes 
of the standing wave, one has more precisely kl = kz + <fi + tt/2. Eq. EH corresponds then 
to the results obtained in ref. Ilol, up to a factor 3/2 due to the fact that the vectorial 

n 

character of the dipole was ignored in ref. [la]. More accurate results for any position of 
the atom and solid angle subtended by the mirror can be obtained from eq. [9""] and 19*71 
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